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ABSTRACT 

We used deep wide-field observations obtained with the Canada-France-Hawaii Telescope to study the blue 
straggler star (BSS) population in the innermost five arcminutes of the remote Galactic globular cluster Palomar 
14. The BSS radial distribution is found to be consistent with that of the normal cluster stars, showing no 
evidence of central segregation. Palomar 14 is the third system in the Galaxy (in addition to w Centauri and 
NGC 2419) showing a population of BSS not centrally segregated. This is the most direct evidence that in 
Palomar 14 two-body relaxation has not fully established energy equipartition yet, even in the central regions 
(in agreement with the estimated half-mass relaxation time, which is significantly larger than the cluster age). 
These observational facts have important implications for the interpretation of the shape of the mass function 
and the existence of the tidal tails recently discovered in this cluster. 
Subject headings: globular clusters: general — globular clusters: individual(PAL14) 



1. INTRODUCTION 

Blue straggler stars (BSS) are hydrogen-burning stars lo- 
cated bluer and brighter than the main sequence (MS) turn- 
off (TO) point in the optical color-magnitude diagram (CMD) 
of star clusters. There is a general consensus in consider- 
ing BSS as the most massive (with M B ss = 1 - 1.4M Q ; see 
Shara et al. 1997, Ferraro et al 2006a) luminous objects in the 
CMD of a globular cluster (GC). Two physical mechanisms 
(both affecting and affected by the dynamical processes oc- 
curring in the cluster) have been proposed for their formation: 
direct st ellar collisions and mass transfer a ctivity in binary 
systems dHills & Davlll976t lMcCreal [J964). The former is 
expected to be particularly important in high-density environ- 
ments, while it should be less efficient in loose GCs and in 
the cluster external regions, with respect to the undisturbed 
evoluti on of primord ial binaries. On the other hand, the dis- 
covery dFerraro et al . 2009) of two distinct and parallel BSS 
sequences in M30, possibly populated by objects generated by 
the two formation processes during the core collapse phase, 
further supports the idea that BSS of different origins can co- 
exist within the same stellar system. Observational proof of 
the connection between the binary and the BSS populations 
in low-density environments is testified by the correlation be- 
tween the number fraction of these two species recently found 
in the core of 13 loose GCs (Sollima et al. 2008). In addition, 
being more massive than the average, BSS tend to sink to- 
ward the bottom of the cluster potential well, under the action 
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of energy equipartition. In most of the surveyed GCs (~ 20 to 
date) BSS appear to be strongly concentrated in the core (see, 
e.g., Figure 2 in Ferraro et al. 2003 and Figure 6 in Ferraro 
& Lanzoni 2009), with their fraction, measured with respect 
to ordinary cluster stars, decreasing at increasing distance 
from the center and sh owing an upt urn in the external regions 
(see also iFerraro eTaTJ [19931 12001 iDalessandro et al.l l2009). 
This feat ure has been interpre t ed by means of num erical sim- 
ulations dMapelli et all 12004 120061: lLanzoni et ail I2007dlbl) 
showing that the central peak is due both to BSS formed in 
place because of stellar collisions, and to mass transfer BSS 
sunk to the centre under the effect of mass segregation. In 
contrast, the rising branch in the cluster outskirts is due to BSS 
generated by the unperturbed evolution of primordial binaries, 
which are preferentially orbiting in regions where the dynam- 
ical friction timescale is longer than the cluster ag e. The only 
two e xceptions currently kno wn are o;Cen tauri ( Ferra ro et alj 
120061) and NGC 2419 (Dales sandroet alj|2008h . where BSS 
show the same radial distribution as that of the other clus- 
ter stars. This fact indicates that these two GCs have not yet 
reached a status of energy equipartition, and their BSS result 
from the evolution of binary systems whose radial distribution 
has not been altered by the the process of mass segregation. 

Based on these results, here we use the BSS radial distri- 
bution to investigate the dynamical state of the remote GC 
Palomar 14 (hereafter Pal 14), located in the outer Halo of the 
Milky Way, at a distance of — 71 kpc (ISollimaet al J 120 111 
hereafter SI 1). 

Pal 14 has been indicated as one o f the best candidates to 
test alternative theories of gravity (Baum gardtet al.l 120051 : 
ISollim a & Nipotjl 1201 Oh and recent analyses have shown a 
peculiar structure and kine matics of t h is clu ster. Adopt- 
ing a sample of 17 stars iJordi et alJ (120091) measured a 
very small veloc i ty dis persion (a v = 0.38 ±0.12 km/s). 
Kiipper& Kroupa (2010) claimed that such a small value is 
not compatible with the presence of the fraction of bina- 
ries (> 20%) expected in a loose GC. They argued that ei- 
ther this cluster hosts a low fraction of binaries, or it consti- 
tutes a "deep-freeze" with an unusually low velocity disper- 
sion. By measuring the cluster mass function between 0.53 
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FIG. 1. — Map of the MegaCam dataset of Palomar 14. The region consid- 
ered in the present analysis (r < 300" ~ 2 rj ~ 8 r c ) is marked with the thick 
solid circle. The two regions contained within the thick solid lines have been 
used to define the sample of Galactic field stars. The tidal radius (r t = 20') 
is shown for comparison and marked with the dashed circle. The iso-density 
contours of the large-scale surface density around Pal 14, showing the tidal 
tails discovered by SI 1, are plotted with thin solid lines. 



and O.78M , Uordi et al.1 (120091) found a significantly flatter 
slope than the canonical value and concluded that either Pal 14 
formed with only few low-mass stars, or it is mass segre- 
gated and lost most of its low-mass stars through interaction 
with the Galactic tidal field. Indeed, two well defined and 
extended tidal tails associated with Pall4, likely due to an 
active process of tidal stripping, have been recently detected 
(SI 1). However, no measurement of the degree of mass segre- 
gation was available to date and the estimated half-mass relax- 
ation time (t r h ~ 20 Gyr; SI 1) is much longer than the cluster 
age (t ~ 10.5 Gyr; Potte r et al]l2008l) . On the basis of these 
results and using N-body simulations, Zonoozietal. (2011) 
concluded that either a primordial mass segregation or a non- 
canonical initial mass function must have been established in 
this cluster after the initial gas expulsion. 

With the aim of clarifying the dynamical state of Pal 14, here 
we study the radial distribution of its BSS population. In Sect. 
2 the adopted photometric dataset is presented. In Section 3 
and 4 the BSS and the reference populations are defined and 
their radial distribution is derived and analyzed. In Section 5 
we present our conclusions. 

2. CATALOGUE AND COMPLETENESS LEVEL 

This work is based on observations performed with the 
wide-field camera MegaCam, mounted at the Canada-France- 
Hawaii Telescope (CFHT). The photometric dataset and the 
adopted reduction procedures are described in S 1 1 . The final 
catalogue includes 40,000 objects sampled in the g' and r' fil- 
ters over an area of 1 deg 2 centered on the cluster. A map of 
the MegaCam dataset is shown in Fig. Q]and the CMD for the 
innermost 300" of the cluster is shown in the left panel of Fig. 

El 

The loose structure of Pal 14 guarantees that the photomet- 



ric analysis is accurate and complete even in the very central 
regions of the cluster. However, in order to quantitatively es- 
timate the level of completeness of the photometric sample in 
the central region of the cluster, we performed a detailed com- 
parison with high-resolution images of the cluster obtained 
with the Hubble Space Telescope (HST). We retrieved a set 
of deep, multi-band images (with total integration times of 
8540s and of 10320s in the F555W and FSUW filers, re- 
spectively) secured with the Wide Field Planetary Camera 2 
(WFPC2, GO-6512; PI: Hesser) We analyzed these images 
by using DAOPHOTII (Stetson 1987). Briefly, an accurate 
Point Spread Function (PSF) was estimated on each frame 
and adopted for the first run of the PSF fitting procedure. A 
master list of stars was extracted from a deep, high signal-to- 
noise image obtained from the montage of the entire dataset. 
Then, the average of the magnitudes measured (through ALL- 
FRAME; Stetson 1987) in each single frame for every master 
list object was adopted as the star magnitude in the final cata- 
logue, and the error of the mean was assumed to be the associ- 
ated photometric error. The WFPC2 catalogue includes 2766 
stars with a mag < 0.2 and sharpness parameter \sh\ < 0.2 down 
to F555W ~ 28, i.e. 5 mag below the MS-TO. The number of 
sampl ed stars is fully compatible with the 2752 stars sampled 
by iJordiet al. (2009) and the completeness study shown in 
their Fig. 4, safely assumed as representative of the quality of 
our WFPC2 photometry, indicates that our WFPC2 catalogue 
is 100% complete well below the MS-TO. The completeness 
of the MegaCam catalogue is quantified as the fraction of stars 
in common with the WFPC2 catalogue in a given magnitude 
interval, with respect to the total number of stars sampled by 
the WFPC2 in the same range. We find 93%, 95% and 100% 
of stars at F555W ~ g' < 23.5,23 and 22.5, respectively. 

2.1. Center of Gravity 

By f ollowing the procedure described in Montegriff o et al.l 
(119951 see also Dalessandro et al. 2008), we have esti- 
mated the center of gravity (C grav ) of Pall4 as the barycen- 
ter of the resolved stars. To this aim, we first performed a 
rough selection of the cluster stars along the canonical evolu- 
tionary sequences in the CMD. Then, through a sigma clip- 
ping procedure, we averaged the a and S positions of all 
the stars contained within three circular area s of rad i us r < 
60", 70" and 80" around the centre quoted bv lHarrisi dl996b . 
Three barycenters were measured in each area by using stars 
brighter than g'= 23.5,23.7,24. The mean of these nine mea- 
sures of the barycenters turns out to be 072000 = 16 h l l m 0.8 s , 
£/2000 = 14°57'27.9", with standard deviations 3.5" and 1.3" 
in RA and Dec, respectively (a large scatter is expected be- 
cause of the extremely low stellar density even in the core 
region of the cluster). This position of the cluster center is in 
agreement (within the errors Act ~ 3.19", A<5 ~ -0.05) with 
the one found by Hilker (200B- 

3. POPULATION SELECTION AND BSS RADIAL DISTRIBUTION 

In order to investigate the cluster dynamical state, we stud- 
ied the BSS radial distribution and compared it to that of 
red giant branch (RGB) and horizontal branch (HB) stars, 
taken as representative of the "normal" cluster population (see 
Ferraro & Lanzoni 2009, for a recent review). This requires 
a proper selection of the samples in radial annuli around the 
cluster centre and an accurate analysis of the contribution due 
to field stars located in the foreground/background of the clus- 
ter. Given the tidal distortion recently detected in Pal 14 (S 1 1), 
its stellar distribution can be considered to have a spherical 
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FIG. 2.— Left panel: CMD of the inner 300" of Pall4, obtained from the 
MegaCam dataset. The selection boxes used to define the BSS, RGB and HB 
samples are marked with solid lines. The average photometric magnitudes 
and color errors are shown by the error bars on the ieit.Right panel: CMD 
of the field region (~ 445 arcmin 2 ; see Fig.l) with the population selection 
boxes superimposed. 

symmetry only for r <, 300" (see Fig. Q] and Fig. 4 in SI 1). 
This corresponds to ~ 2r;, (or ~ 8r c ), a distance where any 
signature of mass segregation in the radial distribution of the 
cluster populations is expected to be well visible. Moreover, 
at these radial distances the field contamination is still accept- 
able with respect to the number of cluster stars. Hence, we 
limit the following analysis to this portion (r < 300") of the 
cluster. 

In order to define the BSS sample we have considere d the 
selection box shown in Fig. [2] While Sandquist (2005|) also 
included stars much closer to the MS-TO (see his Fig. 1), we 
chose a more conservative limit in order to minimize the risk 
of contamination by spurious (blended) sources, which are ex- 
pected to be very few or zero, but could be critical because of 
the small number of stars in Pal 14. The same magnitude range 
was used to define the RGB sample, which include some sub- 
giant stars. The selection boxes used for the RGB and HB 
populations are shown in Fig. |2] We count a total of 24 BSS, 
191 RGB and 24 HB starsQ 

The cumulative radial distributions of the three samples are 
shown in Fig. [3] No significant difference is found among the 
three distributions, thus indicating that none is distinctly seg- 
regated towards the cluster center with respect to the others. 
A Kolmogorov-Smirnov test indicates a 99% probability that 
the three samples are extracted from the same parent popula- 
tion. 

We also used the Anderson-Darling test, which is more sen- 
sitive to the tails of the empirical cumulative distribution func- 
tion than Kolmogorov-Smirnov, to assess possible differences 
at small and large radial distances. Using the ^-sample vari- 

8 We have carefully checked that slightly different assumptions about the 
selection boxes (e.g. using the definitions provided by Sandquist 2005) do 
not change the results of the analysis presented in this paper. 



TABLE 1 
Number of BSS, RGB and HB stars in three concentric annuli 

USED IN THE ANALYSIS. THE ESTIMATED NUMBER OF CONTAMINATING 

FIELD STARS FOR EACH SAMPLE IS REPORTED IN PARENTHESIS. THE 

MEAN RADIUS NORMALIZED TO THE CLUSTER CORE RADIUS IS ALSO 

REPORTED. 
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64(0.12) 9(0.04) 
62(0.47) 6(0.16) 
65(3.11) 9(1.04) 



ant of the Anderson-Darling test dScholz & S tephens 1987), 
we found that for any pair of samples (k=2) the probability 
that they arise from the same underlying distribution is 63%, 
while for the combined set (k=3) of samples, the probability 
raises to 79%. There is therefore no statistically significant 
radial segregation between the samples of BSS, HB and RGB 
stars. 

To further investigate this issue, we studied the radial distri- 
bution of the po pulation ratios by fo l lowin g the procedure de- 
scribed, e.g., in Ferr aro & Lanzonil (120091) . The field of view 
was divided in three concentric annuli, each one sampling ap- 
proximately the same luminosity fraction. The number of ob- 
jects belonging to the three populations was then counted in 
each annular area (see Table[T|i. 

Inorder to estimate the contamination from field stars to the 
selected populations we have taken advantage of the wide ra- 
dial coverage of the MegaCam catalogue. Two areas between 
23' and 25' (i.e., at 3' from r r ) and orthogonal to the direction 
of the tidal tail (thus to minimize the possibility that genuine 
cluster stars fall in field sample; see Fig. Q]) have been selected 
as representative of the Galactic field population. As shown 
in Fig. |2] (right panel) no signature of the cluster stellar pop- 
ulations is found in the field CMD. We count, respectively, 8, 
21 and 7 field stars within the same BSS, RGB and HB selec- 
tion boxes discussed above. We estimate a contamination of 
1.8 x 10" 2 , 4.7 x 10" 2 and 1.6 x 10" 2 field stars per square ar- 
cminute, to the BSS, RGB and HB populations, respectively. 
By taking into account the area sampled by the three radial 
bins, the number of contaminating field stars in each annulus 
and for each population has been derived and it is quoted in 
parenthesis in Table [T] Finally, to compute the average ratio 
of the number of RGB, HB and B SS stars we used Eq. (26) 
in iCervino & Valls-Gabaudl (12003) for uncorrelated Poisson 
variables. 

As shown in Fig. H](upper panel), the radial distribution of 
the ratio between HB and RGB stars is flat, in agreement with 
what expected for normal cluster stars, which all have (essen- 
tially) the same mass. Consistently with what found above, 
and despite their higher mass, also BSS are distributed as the 
normal cluster population and show no evidence of central 
segregation (Fig. HI second and third panels). 

We have also investigated the radial behavior of the dou- 
ble normalized ratio (Fe rraro et alJll993 ). which reports, per 
each considered annulus, the fraction of stars observed in 
a given evolutionary stage (BSS, HB, RGB) divided by the 
faction of cluster light sampled in that annulus by the ob- 
servations: /?p 0p = (A^ s p a r Sp/A^ t p ° p )/(L sam p/Ltot), with pop=BSS, 
HB, RGB. Since the number of stars in any post-MS stage 
scales li nearly with the total lumin osity of the stellar pop- 
ulation dRenzini & Fusi Peccil Il988h . this ratio is predicted 
to be equal to 1 for any (not segregated) reference popula- 
tions, while a larger (smaller) value is expected for popula- 
tions which are more (less) concentrated. In order to estimate 
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FIG. 3. — Cumulative radial distribution of the BSS, RGB and HB popula- 
tions of Pall4 within 300" from the center, with the distance normalized to 
the cluster core radius (r c = 36", from Sll). 

the sampled light, we have integrated the King profile best- 
fitting the surface density distribution (from Sll). The radial 
trend of Rrcb an d ^bss in me same annuli previously defined 
is shown in the lower panel of Fig. @] As apparent, it turns out 
to be constant and equal to 1 not only for the RGB (and the 
HB) stars, but also for the BSS, thus further demonstrating 
that this population has a radial distribution fully consistent 
with that of the reference ones. 

4. DISCUSSION 

The radial distribution of BSS in Pall4 is indistinguishable 
from that of its normal (and less massive) stars. As in the 
case of the only two other clusters showing the same feature 
(wCen and NGC2419: TFerraro et al . 2006; Dalessand ro et alJ 
2008, respectively), this is an observational proof that Pal 14 
is dynamically young, still far from having established energy 
equipartition even in its innermost regions. This is in agree- 
ment with the extremely long half-mass radius relaxation time 
(~ 20 Gyr) recently estimated by SI 1. Moreover, our results 
suggest th at the unusually flat mass function measured by 
Uordi et aT] t2009) cannot be explained by energy equipartition 
developed during the cluster d ynamical evolution, b ut should 
be primordial (as suggested bv lZonoozTe t al. 201 1). We note 
that since the mass range covered by that study is quite lim- 
ited, further investigation is needed. Finally, the negligible 
degree of relaxation of Pall4 suggests that the observed tidal 
tails should not be preferentially populated by low mass stars 
evaporated from the cluster. 

The flat BSS radial distribution also suggests that (as ex- 
pected in such a low density environment) stellar collisions 
played a minor role in generating these exotica and affecting 
the binary population. Hence, as in the case of wCentauri 
and NGC 2419, the BSS we are observing likely derive from 
the evolution of primordial binaries and can be used to get a 
rough estimate of the fraction of such a population. As a first 
consideration, we note that the number of BSS normalized to 
the sampled luminosity (in units of 10 4 Lq, see lFerraro et alJ 



FIG. 4. — From top to bottom: radial distribution of the number fractions of 
HB and RGB stars, BSS and RGB stars, BSS and HB stars, and of the double 
normalized ratios S pop (see text) of BSS (dots) and RGB stars (gray regions, 
with the vertical size corresponding to the error bars). The ratios and their 
errors were calculated assumin g that the number counts of each sample were 
uncorrelated Poisson variables iCervino & Valls-Gabaud 2003). 



11995b is S4 BS s - 2 in wCentauri and S4 B ss ^ 3 in NGC 2419. 
The same ratio in Pall4 rises to 29 (i.e. a value 10 times larger 
than what found in the two clusters with similar BSS radial 
distribution). However this value is not that surprising when 
compared to the field. In fact, as discussed by Ferraro et al. 
(2006), the observed BSS specific frequency Nbss/Nub = 0.1 
in wCentauri turned out to be ~ 40 times lower that what 
observed in the field (Nbss/Nhb = 4) by Preston & Sneden 
(2000). The value found in Pal 14 is Nbss/Nhb ~ 1, in much 
better agreement with the above field sample. 

Under the hypothesis that all the BSS are originated by pri- 
mordial binaries, this possibly suggests that the binary frac- 
tion in Pall4 (and in the field) might be much higher than 
in the other two GC s (in wCen it amounts to /bin ~ 13%; 
Sollima et al. 2007a). A rough estimate of the binary fraction 
in Palomar 14 can be derived from the correlation between 
the measured binary fraction and th e cluster integrated mag- 
nitude found by Sollima et al. (2010) in a sample of 18 open 
and low-density globular clusters. From their Figure 7, and 
adopting My = -4.95 (Sll) we predict / bin ~ 30-40% for 
Palomar 14. In addition, the comparison between the number 
of BSS per unit luminosity (from Ferraro et al. 1995) and the 
fraction of binaries (Solli ma et al.l l2007b) measured in a sam- 
ple of low-density GCs (in which the collisional channel of 
BSS formation is expected to be negligible) shows that while 
BSS-poor GCs (with 8 < 54 B ss < 13) host a small fraction 
of binaries (9 < /bi„ < 16%), the only cluster (Palomar 12) 
with a value of 54bss similar to Pal 14 has a binary fraction 
of ~ 40%. Despite the uncertainty affecting these estimates, 
such a result and, even more robustly, the existence of a non- 
collisional BSS population, suggest that Pall4 might have a 
non-negligible fraction (of the order of ~ 30-40%) of bina- 
ries. 

Deep and high-quality imaging of Pal 14 main sequence are 
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urged for a direct estimate of the binary fraction and a mea- 
surement of the mass function in a wide range of masses. Ad- 
ditional spectroscopic campaigns able to more precisely esti- 
mate the cluster velocity dispersion would also be very useful 
to constrain more tightly the dynamical state of the cluster. 
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